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Abstract 

We study the interpolation from occupation number Fock states to Schrodinger cat states on 
systems modeled by two-mode Bose-Hubbard Hamiltonian, like, for instance, bosons in a double 
well or superconducting Cooper pair boxes. In the repulsive interaction regime, by a simplified 
single particle description, we calculate analytically energy, number fluctuations, stability under 
■ coupling to a heat bath, entanglement entropy and Fisher information, all in terms of hypergeo- 



metric polynomials of the single particle overlap parameter. Our approach allows us to find how 



Ch ■ those quantities scale with the number of bosons. In the attractive interaction regime we calculate 

the same physical quantities in terms of the imbalance parameter, and find that the spontaneous 



symmetry breaking, occurring at interaction Uc, predicted by a semiclassical approximation, is 
valid only in the limit of infinite number of bosons. For a large but finite number, we determine 
a characteristic strength of interaction, U*, which can be promoted as the crossover point from 
coherent to incoherent regimes and can be identified as the threshold of fragility of the cat state. 
Moreover, we find that the Fisher information is always in direct ratio to the variance of on-site 
number of bosons, for both positive and negative interactions. We finally show that the entangle- 
ment entropy is maximum close to U* and exceeds its coherent value within the whole range of 
interaction between 2Uc and zero. 
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I. INTRODUCTION 



The two-mode Bose-Hubbard model is commonly used to describe several systems like, 
br examples, the bosonic double well traps or the superconducting Cooper pair boxes 
3, S, [i1t13|. It has been recently used as the building block for studying spinful bosonic 
systems on double- well lattices [ij]. In spite of its simplicity, being an interacting problem, 
although integrable 15| and whose low energy spectrum was explored in the weak interaction 
regime 16[, it is far from a simple solution for a generic number of bosons and interaction 
strength. Because of that reason, usually one resorts to numerical approaches which are 
simple and very efficient, giving up trying to describe the system analytically, although 
not exactly. A first issue worth being addressed is, therefore, that of providing a simple 
description which allows one to find handy and analytical expressions for several physical 
quantities and their scalings with the number of bosons. In order to discriminate between 
the so-called quantum phase model expressed by occupation number states, and the 
mean field model [l^, formulated in terms of coherent states, emerging from the solution 
of the classical Gross-Pitaevski equation - both models invoked to describe, for instance, 
the charge oscillations in the superconducting Cooper pair boxes - it can be useful to study 
certain stability properties The two models, in fact, behave quite differently in the 

presence of noise induced by a weakly coupled external environment, therefore it is crucial 
to predict how the relaxation time scales with the number of bosons. 

In this paper we propose, therefore, to tackle the problem of describing the ground 
state of a two-mode Bose-Hubbard model from an intuitive and physically transparent 
description, although approximate. Our proposal is based on the observation that the 
norm of a Fock state, which is the ground state for strongly repulsive interacting bosons. 



if 



|0) , (where a}^^ are creation operators on left, L, or right, i?, site) is 



given by the permanent (Per) of a diagonal block overlap matrix 

(?/'f|^f) = Per 





Onxfc 


Ofcxn 


Ifcxfc 



where, schematically, the blocks denoted by "1" are made of all elements equal to one 
and the blocks denoted by "0" are made of all null elements. On the other hand, the 
norm of a coherent-like state, or sometimes called phase state, ground state for free bosons, 
IV'c) = {.iLCi\ + 'Ci?'^ij)"'^'^|0) , with l^ip + l^/jp = 1, is given by the permanent of a fully 
constant matrix 

(V^cl^c) = Per 



Ifixn 




Ifcxn 


Ifcxfc 



A natural expectation is, therefore, that, going from number Fock state to fully delocalized 
coherent-like state, the bosons, initially localized on the two sites, start to overlap, making 
finite the off-diagonal blocks of the overlap matrix. Our ansatz is that the intermediate 
interaction regime can be fairly described by the state I?/'), written in Eq. ([9]), that describes 
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two condensates which can be locahzed on each site in the Fock hmit while merging together 
in the coherent one, and whose norm is given by 

{ip\ip) = Per 





^nxk 


^kxn 





where u is indeed the overlap of the single particle wavefunctions, as we will be seeing. In the 
site-symmetric case, then, this description allows us to interpolate from Fock to coherent- 
like states by varying a single parameter, a;, which has a clear physical meaning and can be 
fixed variationally in terms of the number of bosons and the microscopic parameters of the 
Hubbard model. 

Another important issue is related to the spontaneous symmetry breaking between the 
two wells, which is supposed to occur in the attractive interaction regime. The experimental 
realization of a double-well potential which confines ultracold alkali-metal atoms, is the 



atomic analog of the superconducting Joshepson junction |l7H20|. as predicted several years 
ago [21']. The Josephson equations are valid in the weak interaction regime and predicts 
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25^. 



self-trapping and symmetry breaking of the atomic population in the two wells 
The symmetry breaking can be derived starting from quasiclassical coherent states. We 
will consider, instead, coherent-like states, i.e. I'^/'c); given by the contributions to the 
full coherent states with fixed number of particles, which is the case in real experimental 
situations. Also these states exhibit a symmetry breaking, however we show that a symmetric 



linear combination of two of such states, namely a macroscop.c Schrodinger-cat state 0. 
27] . which has null population imbalance, is energetically favoured, for equal local energies. 
As a result, in the case of fixed number of particles, the symmetry breaking should not occur 
and the imbalance parameter, which does not correspond to the population imbalance, is 
inite for any attractive interaction. This result is in agreement also with numerical results 



281], and, at some extent, consistent with the semiclassical approach in the limit — ?■ oo. 
We find, however, a discontinuity in the imbalance parameter at some critical strength of 
interaction. This is a spurious effect of our ansatz for the ground state, however it reveals the 
fragility of the cat state below a certain value of the interaction. Nevertheless, by numerical 
checks, we observed that, around that value of interaction, some physical quantities have 
large derivatives and the energy changes its behavior. At that point the coherence visibility 
drops while the cat state becomes extremely fragile. The relaxation time, in fact, goes like 
1/A^^ upon an asymmetric coupling with an external bath. One can identify such a threshold 
as the critical interaction below which an infinitesimally small mismatch between the two 
on-site energies can produce a macroscopic population imbalance. 

Finally we investigate the behavior of the entanglement entropy and the Fisher infor- 
mation. The quantum entanglement applied to many body systems has attracted a lot of 



theoretical interest in recent years (see, for instance, Refs. 29|, |30|] and references therein). 



One would expect that the entanglement entropy reaches its maximum in the presence of 



large coherence. We show, instead, that, in agreement with numerical results 28|, the max- 



3 



imum entropy occurs close to the putative symmetry breaking point. In particular we found 
that the entropy exceeds its value obtained in the case of free bosons, in a whole range of 
attractive interaction which goes from zero to twice the critical interaction of the symmetry 
breaking. We finally show that, for almost all the range of interaction, i.e. — oo < U <^ N , 
the Fisher information F is always directly proportional to the variance of on-site number 
of bosons a, as described by the equation F = ia/N"^. 

The paper is organized as follows. In Sec. II we introduce some useful definitions; in Sec. 
Ill we present the model and our ansatz for the approximate ground state in the repulsive 
interaction regime; in Sec. IV we perform detailed calculations for generic two and four-point 
correlation functions, useful to calculate several quantities of interest; in Sec. V and VI we 
recover the exact results, respectively, for strongly repulsive interacting and free bosons; in 
Sec. VII we derive asymptotic behaviors in the large limit for energy, decay rate, number 
fluctuations and coherent visibility; in Sec. VIII we calculate, by our ansatz for repulsive 
interaction ground state, the entanglement entropy and the Fisher information. Section IX, 
instead, is devoted to the analysis of attractive interaction, calculating all the quantities of 
interest as Schrodinger cat state mean values and finally in Sec. X we derive entanglement 
entropy and Fisher information, always in the attractive regime. We summarize the main 
results showing some plots in Sec. XI, and drawing some conclusions in the final section. 

II. GENERAL REMARKS 

Let us first consider N bosons on a lattice with A^^ sites and take the following, not 
normalized, many-body wavefunctions 

N / Ns \ 



where ^ai and rjai are single particle wavefunctions. It has been shown that such states 
are rich enough to exhibit both superfiuid and insulating behaviors [sT]. After defining the 
N X N matrix 




(1) 



(2) 




(3) 



the (A^ + 1) X (iV + 1) matrix 




(4) 
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and the (A^ + 2) x (A^ + 2) matrix 



/ ^ fji fji 

i] ^ii 

\ iln Sim Sml J 



(5) 



where fji is an A^- vector with components rjai, we get the following relations 31 1 



ai ai a] al^\(p) 



(6) 
(7) 
(8) 



where Per(y4) is the permanent of the matrix A. Quite in general, since the calculus of 
permanents is not easy, it is hard to find handy and analytical expressions for the correlation 
functions. 



III. THE INTERPOLATING STATE 



Now let us take N = n + k bosons on two sites, like in a double well, and consider the 
following many-body state 



Gl^+^RoIj) (^^2La{ + ^2Ra^R^ |0) 



(9) 



which is a simplified version of the permanent in Eq. ([T]), where Va such that a < n, S,ai = 
while Va such that n < a < {n + k) , = ^2i- Defining the unit n x k matrix (with n rows 
and k columns) with all elements equal to 1, 



j(n,fc) 

we can construct the overlap matrix 



(\ 1 ... l\ 



11... 1 



(10) 



where 



&= |6L|' + |eiij|'= |6L|' + |e2ij|' 



(11) 



(12) 



is the normalization of the single particle wavefunction, therefore we can set 6 = 1. However, 
for the moment, in order to derive more general relations, we will keep the writing h. The 
other parameter is 

^ = i2LilL^i2RilR (13) 



5 



which is, therefore, the single particle overlap term. For normalized single particle wave- 
functions, namely when 6 = 1, we have < |tj| < 1. 

Eq. is the simplest permanent state which interpolates, keeping the same functional 
form, from number Fock state (when, for instance, ^il = ^2R = 1 and ^2L = = 0, such 
that we have n bosons on the left well and k bosons on the right well) to a coherent-like 
state (when = ^2L and = ^2r)- 

The four parameters, C,il,^2L,(,ir,^2R, are not independent, but linked by the two nor- 
malization conditions, Eq. (fT2l) . On the contrary Eq. (fT3|l is not a further constraint but 
leads to a different parametrization, where the overlap coefficient u replaces one of the ^- 
parameters. Moreover, from Eq. one can consider n as another free parameter while k 
is fixed hj k = N — n. Therefore, at the end, we have three free state parameters (in general 
two complex numbers and one integer): ^i^, (one among the four ,^'s), the overlap u and n. 

These values can be related to the parameters of the microscopic model. Let us introduce 
the Bose-Hubbard model 

H = (^a{aR + a^Oi j + n{nL- ur) + ^ (^^(nL - 1) + nR{nji - 1)) , (14) 

where t is the hopping parameter, /i is the difference of local chemical potentials induced for 
instance by an mismatch between the two wells, and U is the on site interaction. One can 
therefore link the three parameters .^i^, u, n with the three parameters of the Bose-Hubbard 
model, t, /i, f/, by variational analysis. Actually, defining 

E[ei,,n,c.|t,/.,f/] = i||M, (15) 

one can partially fix the set of state parameters, requiring that S^^^^E = 6^E = 6nE = 0, 
establishing, in this way, a connection between the set of the state parameters with the set 
of the Hamiltonian ones. The residual arbitrariness can be removed if we require that the 
amplitudes are real. 



IV. CORRELATION FUNCTIONS 



In order to calculate two and four points correlation functions we have to determine 
the permanents of the matrices Q, D and /. As we will see, all these permanents can 
be calculated analytically. Let us start considering the normalization of the many-body 
wavefunction Eq. ([9]). This quantity is given by 

{^P\^|J) = Per(fi"'^) = n\k\ 6^^+^ 2^1 (^-n, -k, 1; (I6) 

which is derived by induction, applying the Laplace theorem for the permanent and where 
2F1 {a,b, 1; z) is an hypergeometric function. Since the first arguments are integer numbers. 
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2-Fi is a polynomial, particularly, using the definition of the Jacobi polynomials P!i''^\z), it 
can be written as follows 



2F1 (-n, -k, 1; z) = pio^-"-k-i)(^i _ 2z) . 
Actually, one can verify that Per(fi"''^) satisfies the following recursive equation 



Per(fi"''=) = n6Per(fi"-^''^) + k\ \uj\'b' 



|2j,fc-2 



k-1 



n 



! 6" + n ^ -— Per(r]"-i'^ 



(17) 



:i8) 



and check that Eq. f[T6|) is a solution. In order to calculate two and four point correlation 
functions it is convenient to calculate the permanents of the following matrices 



n,k 
lu 



Jjj{n-l,n) ^* j{n-l,k~l) 
^J{k,n) ljj{k,k^l) 



o 



n,k 
Id 



n-1) 



■ j(n,k) 



^j(k-l,n-l) ^j{k-l,k) 

with dimensions {n + k — 1) x [n + k — 1) and 

^j{n-2,n) ^*j{n-2,k-2) 
^J{k,n) J^j{k,k-2) 



(19) 
(20) 



n,k 
2u 



o: 



n,k 
2d 



})j{n,n-2) j(n,k) 
^j{k-2,n-2) ^j(fc-2,fc) 



(21) 



(22) 



with dimensions {n + k — 2) x {n + k — 2). Their permanents, together with Eq. (fT6ll . 
are the building blocks useful to construct the permanents of the D-matrix, for the two- 
point correlation functions, and of the /-matrix, for the four-point correlation functions. By 
recursion, after several algebraic steps, we get 

k-l 



Per(C>;^*=) =n!(A;-l)!a;6"+'=-2^ ^fA\ 

(i—n V 



n, 1; 



£=0 



^er{Ol^) = {n- l)\k\ tu* 6"+'=-' ^ 2F1 ( 1 - k, 1; 

n-2 

k{n -2)\{k- 2)1 

m=0 

fe-2 

n(n-2)!(A;-2)!w*2 5"+^-4^ 



\IjJ\ 



Per(a 



n,k\ 
2d ) 



1=0 



1 + m) ^ 2^1 ( - 

n-2 , 



■m, — 1; 



1^ 



■m, 



' U2 



However we will need only a couple of those permanents since 

Per(0;^^) = Per(0;^^'=)* , 



Per(02";') = Per(02T) 



n,k\ 



(23) 
(24) 

,(25) 
.(26) 



(27) 
(28) 
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We now write the [n + k + 1) x [n + k + l)-matrix 



D, 



n,k 



(29) 



where i,j = L,R and with = ^^^(1, 1, 1), a ra-vector and = ,^2^(1) 1); a fc- 
vector. Its permanent, related to two-point correlation functions, can be written as follows 



|aia]| 



PeiiD'^j") = 5,,Per(n"''=) + n^u f<,Per(fi"-i-'=) + kC2,Per{OZ') 



+ki 



2i 



kC2iPeT{Q'''''-^) + <*,Per(0 



n,k\ 
'Id ) 



(30) 



Now defining the [n + k + 2) x [n + k + 2) matrix 



jn,k 
lijm 



^J{k,n) ^j(k,k) 



(31) 



with I, i,j, m = L, R we get 
{ip\aiaia\a. 



6i^Per{D^f) + 5,,Per(L>rr ) + S^mPeT{D 



ij J - -ij 

n— l,fc— 1 



5ij5imPei{U 



n,k\ 



+{k - 1)6. ( (n - mMO^;,') + kC^MO 



Id 



+ kC^m^eriD^J 



n,k—l^ 



+<1Z 



5 2m 



5..Per(0^;') + k^2^ (<,Per(fi"-i'^-i) + {k - m.^eriO^,;,'-')) 



+ in - l)^u { (k - l)aPer(02";') + <,Per(Or. 



5 2j 



<LPer(D: 



n— l,fe\ 



.(32) 



Now we are in the position to calculate analytically, in terms of the Jacobi polynomials, 
several quantities, like the total energy, the number fluctuations, the decay rate when the 
system is coupled to an external environment and the coherence visibility. 



A. Energy 

Supposing that our system is described by the Hamiltonian given in Eq. f |T4l) . we can 
calculate the energy of our many body state, Eq. ([9]), as follows 



E={H) 



t [PeT{Dli) + PeT{D-/,)\ , f/ ^Per(/2i\j + Per(/ 



I Per(fi"''=) / "2 I Per(n"''^ 



RRRR / 



+ (^ _ 2^) Per(D2f) _ Pei{Df^) 

^Per(n"''=) ^Per(fi"''=^ ^ ' 
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B. Number fluctuations 



In the same manner we can calculate the charge fluctuations, for example on the left site, 

^L = (ni)-K)^ (34) 
as the variance of the number operator, where 

- 7-J—\ - n„„^On ^\ ^ ^ ' 



Analogously one can write the charge fluctuations on the right site, namely ctr. 



^^^2 ^ / ^^14^^ 1^)^ = ( _ 1 I . (36) 



C. Decay rate 

In this subsection we suppose our system weakly coupled to an external environment. 



therefore the total Hamiltonian is 32 1 



Ht = H + Hb + X (wB + W^B^^ , (37) 

where Hb is the environment Hamiltonian while the last term is the weak coupling between 
the bath and the bosons, i.e. A <^ 1. is a bosonic operator. The heat bath acts as a 
source of noise and dissipation. We assume the environment to be in an equilibrium state 
and that behaves like a white noise, {B{t)'^B)E — 5(t). The presence of the bath leads to a 
master equation for the density matrix, p = ^-^^^^ of the Kossakowski-Lindblad form 33 1 



dtp=-i[H + H^^\p\+V[p\, (38) 

where 

V[p] = 7 (wpW^ - ^{W^W, p}^ + 6 (w^pW - ]^{WW\ p}^ 

+/3 [wpW - \{W\ p}^ + /3* [w^pW^ - \{W^\ p}^ , (39) 

with 7 = A2/+^c/x(5t(x)5)£;, b = J^^ dx{B{x)B^)E, and (3 = J^^ dx{B{x)B) e, 
some coefficients, /3 G C, 7 > 0, 5 > 0, satisfying the condition 7^ > which ensures 



the complete positivity (see Ref. [32| for more details). These coefficients incorporate the 
dissipative effects on the dynamics. The term H^'^\ is an environment induced additional 
Hamiltonian term whose explicit expression is not important for our purposes. Here, in fact, 
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we focus our attention to the stability of our initial pure state 
the constant decay rate defined by 



at t = 0, calculating 



where 



{ip\w^w\ij) 



-fT^ + STs + 2Re[/3T 



/3 



{il)\WW^\%l)) 
{il)\WW\%lj) 



\w\ 



{Mw\ij) 



^l)\W\il} 



(40) 

(41) 
(42) 
(43) 



Let us choose as a generic single particle operator, linear combination of an hopping term 
and right and left density operators 



W = CLnL + CRnR + Ch a^j^ai ■ 
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(44) 



I.e. 



Usually only the hopping operator is considered in the coupling with the bath 
cr = cl = 0. Here, instead, we consider a more general operator. Without loss of generality 
and without spoiling the complete positivity, we could put /3 = and consider only or 
to show how the decay rate scales with the number of bosons. We have therefore 

= |cl|Vl + |cr|Vk + Ic/^pr^ + 2Re(cLC^) ((n^n^) - (ni)(n^)) 
+ CLcl{{nLa)^aii) - {nL){a{aR)) + Chcl{{a'^aLnL) - {nL){aRaL)) 
+ CRcl[{niia\aR) - {nR){a\aR)) + ChC*R{{a'^aLnR) - {nR){a'j^aL)) 



(45) 



where 



..) = ^^l^j^^^ , ctl and aR are defined by Eq. and 



{a^L^RaRai) - \ {aLaR)\ 



LL, 



Per(fi 



Per(L) 



n,k\ 
LR) 



Per(^]"''=) 



(46) 



which is the common 7-contribution to the decay rate when the bath induces only hopping of 
particle between the two sites ll| . When the environment couples also to the local densities 
we have to calculate, in addition. 



{ninR) - {nL){nR) 



P^'^i^LRLR) 



FeT{Dlt)PeriDl 



RRJ 



Per(fi 



{nLa[aR) - {ul) {a[aR) 



{nRa{aR) - {nR){a{aR) 



Per(/ 



n,k \ 

lrll) 



Per(D"''^ 



LR) 



n,k\ 



Per(n 
Per(/^^iK) 



Per(fi 



n,k'\ 



yPer(fi"''=) 



- 1 



Per(D 



rl) 



n,k'] 



- 2 



PeriD'H) ( Per(D^i 
I Per(fi"''=; 



Per(fi"''=) Per(fi"''=) 
and conjugate terms. Analogous calculations can be done for and 



1 



Ferial 



n,k\ 



(47) 

(48) 

,(49) 
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D. Visibility 



In cold atom physics, in order to detect the coherence properties of the condensates, 
the quantity which are commonly used is the momentum distribution, namely the Fourier 
transform of the one-body density matrix C{x,x') = (a(xVa(x')) , where the brackets mean 



the ground state average. It has been shown [34|-|36| that the momentum distribution, 
n{p) = J dxdx'e~^^^^~^'^C{x, x'), can be written as n{p) = no(p)(l + acos{pd)), where no{p) 
is the momentum distribution in the incoherent regime which depends on the details of the 
feasible double-well potential, d the distance of the two minima of the double- well and a the 
so-called visibility. In terms of our site-operators, the visibility can be defined by 

2\{a{an)\ 



a 



N 



which, if valuated using our permanent states, is simply given by 



a 



2 

N 



Per(L) 



n,k\ 

rl) 



Per(fi 



(50) 



(51) 



V. FOCK LIMIT, a; = 

Here and in what follows we put b = 1, namely we fix the normalization of the single 
particle wavefunctions. For a; = we have Per^Oiu) = Per(Oid) = PeT{02u) = Per((9d2) = 
while 2F1 (— n, —k, 1; 0) = 1, Vn, k, therefore 



n,k\ 



Per(fi 
Per(/^r/') 



nlkl , 



(52) 
(53) 



+ {k' - k)UC2M*2m + nki^^u^2^ + ^U^2l) [^^^2, + Cl,C2m 

The charge fiuctuations, then, reads 

CTL = n\^,L\'\U\' + m2L\'\^2R\' + 2nk\U\'\U\' 
and the decay rate F^, given by Eq. fHBl) . is simply 



(54) 



(55) 



(56) 



where, of course |CqlP + = 1 is understood. So far we have used only co = 0, which is 

a necessary but not sufficient condition for having spatially separated Fock state, as we will 
see in Sec. VIII. A. 4. Now, choosing the Fock state as follows 



eiLai)"(6i?a5?)'|0) 



(57) 
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we are forced to put = SiL and = in Eqs. (153 ||5^ . This implies that the occupation 
numbers are {ul) = n and (nn) = k, and the energy, the charge fluctuations, the visibility 
and the decay rate are simply given by 

E=j{n{n-l) + k{k-l))+ix{n-k), (58) 

CTL = CTK = , (59) 

a = , (60) 

T-^ = \ch\^nik + l), (61) 

and, analogously, Ts = \ch\'^k{n + 1) and F/j = |c/ip(r;, + l)(fc + 1). The decay rate, in the 
Fock limit, is, therefore, oc N'^/A, for n ~ A; ~ N/2, and is due only to the coupling of the 
bath with the hopping-like term. 



VI. COHERENT LIMIT, lo = 1 

For a; = 1 we have that 



(-n, -fc, 1; 1) = I + ^ 1 , (62) 



n 



2Fi(l-n,-^,l;l)= ( ''^^ ^ ) , (63) 



n-2 k~2 



EN ^ / N n(n — l)(n + k — 2\ 

il + m)J2 2FA-m,-i,l-l) = ^—^[ . (64) 

m=0 i=0 H \ n J 

We get, therefore, the following simple expressions 

Per(^]"'^) = (n + A;)! , (65) 

Per(0;^'^) = Per(0;^^'=) = {n + k - 1)\ , (66) 

FeriO-,;!') = Per(02"^) = {n + k - 2)\ . (67) 

The coherent state is obtained when 

= = 6l , (68) 

U = ^iR = 6ii , (69) 

so that Eq. ([9]) reduces to 

\^c)= {^Lai + ^Ra^^Y \0), (70) 

where = {n + k), the total number of bosons, and + |^_rP = 1. In the coherent case 
we obtain, therefore, simply 

Per(D:;.'^) = N\{6,,+NU;} , (71) 
+{N' - N)UiC,CJi ■ (72) 
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Notice that Eqs. (17T1) and (172|) are valid also for a lattice, namely, in the presence of many 

sites, when the state is given by 1^) = (^f'^i4) |0). From Eq. (^Sj we get that the 
occupation numbers are the following 

K)=iV|aP, (73) 
{nR) = N\U'- (74) 
Energy, charge fluctuations, visibility and decay rate are given by 

E = -l + MI) N + /.iVd^r - + I (lar + NiN - 1) , (75) 

a = aL = aR = N\^L\'\U\ (76) 

a = 2\^RCL\, (77) 

r, = N{\cL- cr\' + \ch\' \U' + 2Re[(ci|eL|' - cr\U')cIM1]} , (78) 

and analogously, Ts = N {\cl - cr\^ + \ch? 1^1' - '^M^ClI^l? - crIU'KUCl]}- 

What we found is, therefore, that, in the coherent regime, both the number fluctuations 
and the decay rate scale linearly with the number of bosons, no matter how the bath is 
coupled with the bosons. The couplings affect only the prefactor of the decay rate. Before 
we conclude this section we present here a couple of interesting effects which are supposed 
to occur in the coherent-like regime. 



A. Self-trapping 

Let us consider small interaction so that the coherent state still well approximates the 
ground state or suppose that we can prepare the system in such a state. After defining the 
variables z, the relative charge imbalance, and 0, the phase difference between the left and 
right amplitudes, 

z^{{nL)-{nR))/N=\^L\'-\U\ (79) 



e 



we can rewrite the total energy Eq. fl75|) . as follows 



(80) 



E = ~NVl-z^ cos(p+ ^N{N - 1)(1 + z^) + ^Nz . (81) 

Eq. (18T|) is the Josephson energy, where the hopping term is proportional to cos0. The 
effective Hamiltonian describing the evolution of z is = E/N which gives the following 



equations of motion [23|, IstI 



^Vl-z^ sin0, (82) 

Uz,^^ . tzcosd) 

— -a. (83) 
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(84) 



The Josephson current is, therefore, Ij = Nz = — z'^ sin 0. For /i — )■ 0, we get the 

following non-trivial fixed points 

^ ~ ^ U{N-1) ' 
= 0,±7r, (85) 

which means that, for a particular imbalance, the particles do not arrange themselves to 
reach the symmetric configuration. 



B. Attractive interaction: Symmetry breaking 



We now consider the attractive interacting case, U < 0, supposing that the state remains 
coherent also for negative U, not only strictly for zero interaction. The total energy of a 
coherent state is given by Eq. (1751) . The energy profile, for fi = and = 0, has only one 
minimum at 

1 i 

(86) 



1 

2' 



laP = ^, for 



while develops two minima at 



2^ 2U{N-1) ' ^<T3^ 



(87) 



For U < Yz^, then, the ground state becomes twice degenerate and the charge imbalance, 
z = 2\C,l\'^ — 1, is finite, given again by Eq. (E 



VII. LARGE N LIMIT, uj G (0, 1] 

In this section we derive the asymptotic behavior of the charge fluctuations and of the 
decay rate in the large N limit. The result is valid as long as 

\uj\ > A^"^ (88) 

and for k = n = N/2. In this case we obtain the following asymptotic behaviors for the 
permanents appearing in the correlators, Eqs. ( I30l l32|) . 

Per(l]"-^'") 



Per(fi"'" 
Per(Ci;") = Per((]"''^) 

Per(02T) = Per(l]"'") 
Per((]"-2''^) = Per(r]'^'" 
Per(fi"'^'"-i) = Per(fi"'" 



1 - 


u 


+ An 


4(1 + 







- 1 + 


4 


u 


n 


4uj*{l + 







+ 0{n-^) 
+ 0{n-^ 



2a;2(|^ 


-1 + 


UJ 


n) 


\liJ 


3(1 + 


uj\Y{2n 




l)ri2 



2 


- 2 




+ 2n 


(1 + 


u\y{2n - 1)^2 



(1 + |w|)2(2n- l)n 



+ 0(n-^) 
+ 0(n-^) 
+ 0{n-* 



(89) 
(90) 
(91) 
(92) 
(93) 
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and Per(C>;^") = Per(C»i;")*, Per(n'^'"-i) = Per(n"-i'"), Per(f]"'"-2) = Per(f^"-2'"). Notice 
that for w = 1 only the leading terms survive and we recover the exact results for the 
coherent states, Eqs. (165116 7p . 

Using these relations we can calculate the charge fluctuation and the decay rate. For 
simplicity we now choose all ^'s real so that a; is also real. Since we are considering the 
simple case with k = n, namely a symmetrically occupied double well, we can use the 
following parametrization 



U = = ^ ^Jl + Vl^ , (94) 



6l = ^iR = 1 - v^l^' (95) 

which fulfils the normalization conditions and the definition of the single particle overlap. 
We have also checked that, for fi = and repulsive interaction, U > 0, the balance con- 
ditions [k = n and = ^2Ri = C,2l) naturally minimize the total energy. After this 
parametrization we can calculate, for instance, a and F^, and, after expanding in A^~^, we 
obtain, at the leading order, 

(jl = (yR = (y ^ — , (96) 



F ~ — 



N - 



and, analogously, F^ ~ 



4 

the asymptotics written a 



|2 



(97) 



^ 'y/uj{cL — cr)| . Calculating the total energy, Eq. fl5^ . with 
Dove, Eqs. flHUjl -f l^ . we obtain 



U(8 + Auj{uj^ -3)-2N + 2uj{A + 00- 2uj^)N + N\l + oofiN + a; - 3)) 

+ ^(l + u;'-2a;(iV + l)) , (98) 
and imposing S^oE = 0, we get the minimum energy when 

c^o ^ -T^= ' (99) 

Vt + NU 

which is consistent with the limit u ~ N^^/"^ ^ A^^^ the asymptotic is based on. As a 
result, the decay rate at f/ = 0, is 



A^ 

T^,:,-\CL-CR + Ch\^ (100) 



and F^ ~ ^ \cl — cr — Chf- For U/t of order one and for large A^, instead, the leading term 
is given by 
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We know that ~ A^^|c/ip/4 for the Fock states, meaning that when lOq ~ A^~^, which 
imphes U ~ tN , namely, when the interaction is of order A^, we enter the Fock regime. 
The charge fluctuations, far from the Fock regime, are 



UJoN 1 ft ^./r, 

4 4V f/ 



(102) 



For ?7 = 0, namely = 1, we recover the resuh of Sec. |Vll with = j^^l = 1/^2, smce, 
without interaction, the coherent state is the ground state. Finally, making an expansion in 
N-^ of the total energy, Eq. ( !98|) . calculated at the minimum, Eq. ( l99ll . we get 



t 



Eo = --N +-UN{N -2) + -VtUN + 0{l), (103) 

which, in spite of its simplicity, is a very good approximation of the ground state energy 
for intermediate interaction, i.e. <^ ^ <^ N . By exactly diagonalizing the Hamiltonian 



for different values of A^, in fact, one can easily check that Eq. fllOSp deviates from the 
exact ground level Bex by an error, Eo — Eex, of order 0(1) ~ (I + (relative error 
of order 0(iV~^)). Notice that the coherent energy, Ec = -|iV + ^N{N - 1), deviates 
from the exact one by an error 0{N) ^ and the Fock energy, Ep = jN(N — 2), by 
an error 0{N) ^ (both relative errors of order 0{N~^)). As a result, Eq. f ll03p can 
be considered as a non-perturbative analytical expression of the ground state energy for 
intermediate interaction. 

Let us finally, consider the visibility, Eq. (I5T1) . This quantity is equal to 1 in the coherent- 
like state {U = 0) and decreases by increasing the repulsive interaction {U > 0). Using Eq. 
(I30D and Eqs. (189]), (190]), for ri = A; = A^/2, and Eqs. (El]), (|95]), we obtain, in fact, for 
< U/t N, at the leading orders, 

a ~ 1 - ^ " ~ 1 - ^ ^ . 104 

This result for the coherence visibility is in a very good agreement with numerical results 



28l |. To conclude, both the energy, Eq. f ll03p . and the visibility, Eq. fll04p . have been 
successfully compared with exact diagonalization results, validating therefore our ansatz of 
the interpolating state in the regime of repulsive interaction. 



VIII. ENTANGLEMENT FOR > 

In this section, we conclude our study 
Fisher information of our permanent state. 



calculating the entanglement entropy and the 
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A. Entropy 



Given the wavefunction ([9]), we can derive the reduced density matrix for the left site, 
for instance, by tracing out the right one (for a review, see Ref. 29| and references therein) 



1 



P 



Tr«(l^)(^l) 



(105) 



where Tr^j is the trace over the right site. More exphcitly, we have to perform the following 
summation 



n+k ^ 



This is a diagonal (n + fc + 1) x (n + fc + l) matrix, 
can be written analytically as follows 



(106) 

Pi5ui, whose diagonal elements 



Pt 



for < £ < /c and 



|e2Lp"|eii?P'l6iJp('-') {n + k- £)! k\ I2F1 (-n, 1 - £ + A;; 



n! £! ((A; - ^F, ( - n, -k, 1; + i2RilR\) 



(107) 



Pi 



\iiR\'^^l\n\ 2Fi(^i-n-k,-k,l + i-k;^^ 



r(2L 



k\{n + k- £)\ {{£ - k)\y 2Fi( - n, -fc, 1; + ^2rCir\) 



(108) 



for k < i < n + k. The ^'s are related by normalization conditions, Eq. fll2p . We can 
therefore calculate the von Neumann entropy 



n+k 

S = -'^Pe log2 , 



(109) 



with log2 the logarithm to base 2. We have shown that also the reduced density matrix and 
the entropy can be written in terms of Jacobi polynomials of oo and of the ratio . 



1. Fock limit 



In the Fock limit, a; = 0, choosing = diR and = ^iL (the other choice is equivalent 
but one should take care of the ratio , the argument of the Jacobi polynomial appearing 
in the numerator of the density matrix), from Eqs. f ll07p . (11081) . we get 



Si 



(110) 



As a consequence the entropy is simply S = 0. 
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2. Coherent limit 



If = ^2L and = ^2R, then w = 1, and the reduced density matrix becomes simply 

n' 



pi 



;iii) 



where, of course, N = n + k and = (1 — In this case we can calculate the 

asymptotic behavior of the von Neumann entropy for ^ 1, since the binomial distribution, 
Eq. fillip , approaches the gaussian one 



Pt 



exp 



{i - iV|6 



|2\2 



;ii2) 



obtaining the following asymptotic behavior for the entropy [31| 



(113) 



3. Intermediate case, with n = k 



Now let us consider the case of = n = N/2 and Eqs. flMj) . f l95|) . Applying Eqs. f ll07p . 
f llOSp and f llOOp we get an entropy which is well approximated, for w ^ 1/A^, by 



(114) 



as one can see from Fig. ([T]). At a; = Uo, for U/t of order one, we get 5 ~ ^ log2(A^)+ const. 




FIG. 1: (Color online) Entanglement entropy, S, in the repulsive regime, as a function of w, for 
N = 300 (upper solid line) and = 50 (lower solid line). The dashed lines are the approximated 
behaviors given by Eq. ()114p . for N = 300 (upper dashed line, which cannot be distinguished from 
the solid line) and = 50 (lower dashed line), which deviates from the solid line only for small uj. 
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4- A special case: an almost maximally entangled state 

In this paragraph we will consider a simple state which has the same functional form of 
Eq. ([H]), with same amplitudes in modulus, obtained by a phase deformation of the balanced 
coherent-like state. We will show that the phase can be tuned in order to get an almost 
maximally entangled state. Let us consider, therefore, Eq. (|9]), with 

eiL = 6l = ^iR = e-'*^2R = 1/ V2 , (115) 

such that the not normalized state reads 

1^) = {a[ + a^^y {a{ + e^^a^^)' |0) , (116) 

then the overlap parameter is given by 

u; = i(l + e^^). (117) 

When = we recover the coherent state. Using Eqs. fll07p . f llQSp and fll09p one can verify 
that the state in Eq. flll6p is more entangled than the fully delocalized coherent state, i.e. 
S{(f) 7^ 0) > S{(f) = 0), as we can see in Fig. [2J 




1 2 3 4 5 6 





FIG. 2: (Color online) Entanglement entropy for the state in Eq. (jll6p . with n = k = N/2, as a 
function of (j). For = 300, the entropy reaches its maximum at = vr it 0.02 which is S{(/)) ~ 8.13. 
The dashed line is the upper limit entropy, log2(A^ + 1), which, for N = 300, is ~ 8.23. Exactly 
at (p = TT there is a narrow dip where S{(f) = vr) ~ 6.96, nicely approximated by Eq. (jl22p . which 
gives w 7.23. At (p = and 2-k, the entropy coincides with Eq. (jll3p . giving S{(j) = 0) 5.16. 

In particular, for = tt the single particle overlap is a; = 0, like the Fock state. Defining 

aI = (4 + 4)/v^, (118) 
4 = (4-4)/v^. (119) 

the state in Eq. (11160 is, indeed, a Fock state in the new representation, {tp) = {o\)"' {(12)^10), 
very unstable under coupling to an external bath, in fact F^: ~ n/c/2 (see Eq. ( 156|) ). and, at 



19 



the same time, with a quite large number fluctuations, cr ~ nA;/2, (see Eq. (155|) ). For = vr 
and k = n = N/2, the reduced density matrix is simply given by 

{2iy.{2n - 2i)\ 
4" (£!)'((n-£)!) 
P2^+i = . (121) 



Studying the proflle of the reduced density matrix, Eq. (11201) . one can check that, it can be 
approximated as p2e ^ 1/n = 2/N, yielding, therefore, the following approximated value for 
the entropy 

^^log2(iV)-l. (122) 

The state Eq. (11161) with (p = tt and n = k, is, therefore, almost double entangled with 
respect to the coherent state (0 = 0), as one can see by comparing Eq. (I122p with Eq. (11131) . 
In conclusion, the state in Eq. (I116p . at ~ vr, is almost maximally entangled since the 
entropy approaches the upper limit for N bosons, i.e. log2(A^+ 1), as shown in Fig. [21 What 
shown is a simple example of how the entanglement can increase by losing coherence. 



B. Fisher information 



The Fisher information is deflned by 28|, |38|-|40| 



F= J^2 ' (123) 

which can be written as follows 

F = ^{<^L + <yR-'2{{nLnR)-{nL){nR))] . (124) 

If we consider the case with n = k = N/2, after choosing ^ real and parametrizing them as 
in Eqs. (El]), ([95]) we get 

F ^ ^ , (125) 

for u <^ N^^, and 

F^^, (126) 

for u ^ A^"^. The overlap of the bosons can be seen, therefore, as the parameter which 
encodes the quantum information. In the latter case, i.e. u ^ A^~^, we flnd, in fact, that, 
at least at the leading term, 

{riLTiR) - {nL){nR) -a (127) 

where a = = or. Equation (I127P becomes exact in the coherent limit, a; = 1. This 
implies a very interesting relation, valid at least for u ^ A^~^, 

(128) 
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For finite interaction and at the ground state energy, we know that Uo — \Jt/ (NU), therefore, 
we get 




(129) 



We notice that the Fisher information and the decay time = l/F^, for finite U/t, behave in 
the same way, namely they are proportional, i.e. F ~ r-y|c/ip/4, at least for a symmetrically 
occupied double well and large N. We can say, therefore, that the greater is the Fisher 
information, the longer the state survives under coupling to an external environment. 



IX. ATTRACTIVE INTERACTION: THE SCHRODINGER CAT 



In this section we will study the evolution from a coherent-like state to a NOON state, 
namely \N)l\0)ji + \0)L\N)ji, reached in the limit U — )■ — oo. As we have seen before, this 
regime seems to be characterized by a symmetry breaking, as soon as t/ < Uc, where 

We will consider, therefore, the following "cat" state, superposition of two unbalanced 
coherent-like states, 

M = M + \^Pr), (131) 

with 

\^L)={eai + ^^a^nY\^)^ (132) 
\i^R)={^^4 + ea^RY\0). (133) 

where |^^| > l^'^l, and + = 1, meaning that is more weighted on the left site 
while {iPr) on the right site. After defining the charge imbalances for lipi) and {iPr), one 
opposite to the other, 

z = ZL = -zn = \e\'-\n\ (134) 

the energy in the quantum coherent regime, can be written as in Eq. ( IHTl) (with /i = 0). 
Imposing dE/dz = and dE/dcf) = 0, we get two minima, for any cj) = 2mT, 

z - ^VUWEIIEI n 35) 

UiN-1) ' ^^^^^ 



as soon as U < Uc- The two solutions correspond to the state or to the state \iPr). 
Actually, the energy is even reduced, as we will show in what follows, by taking the symmetric 
linear combination of those two states, which leads to the cat state written in Eq. 01311) . It 
describes a superposition of a cat sitting on the left site but with his tail on the right one 
and viceversa. When the cat withdraws the tail we get the NOON state. 
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In order to calculate the quantities of interest on the cat state, we first write the following 
correlators, with {ips) = {I^l), I^a)}, 

{^.\^,) = N\bl,^Ferin'''%s,s')), 



(136) 



(137) 



{iJs\aiaia\al^\iJs') = SimPeT{Df-'^{s, s')) + 6ijPeT{D^J^{s, s')) + 5irnPer(L);^'°(s, s')) 



5,/,„Per(fi^'°(., .')) + N'^tOeriDf-'''is, s')) 



where s = L,R and 



R-e 

L _ cR — c< 



^R — ~ S 
&L,L = bR,R = 1 , 

bL,R = bR,L = e^^*+^^e*- 

The normalization of the cat state is, therefore, given by 

= 2N\[i + {e^^* +^^eY]=^N\ i + (^yr3^cos, 

2Nl 1 + ( VT^ 



(138) 

(139) 
(140) 
(141) 
(142) 



AT 



cf>=2mT 



N 



(143) 



Without loss of generality, choosing ^< real, from Eqs. f ll37p . fll38p . we get the following 
quantities, useful to calculate energy. Fisher information, number fluctuations and decay 
rates, all in terms of z, the imbalance parameter. 



@ 



{nL)@ = {nR)m = 



{nl)@ = {n\) 



@ 



{il)@\a\aL\i^@) _ N_ 
{^@\tp@) 2 

(V^@|a[aLa[aL|V'@) 



N 



N / ' 



(^@|a[aLa^a/j|?/'@) _ 



(iV-1) 1 



1+ (v^r^) 



N \ ' 



^ \% O-RIi 

ttRja 



(144) 

(145) 
(146) 

(147) 

(148) 
(149) 



From the equation above we can now write the analytical expression for the energy as a 
function of z 



1 + (v^r^) 



N 



+ ^N{N-1) I 1 + 



1 + (v^r^) 



N 



(150) 
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If z = we recover the energy in the coherent-hke state, namely E = —tN/2 + UN{N — l)/4:, 
while for 2 = 1 we get the energy in the NOON state, E = UN{N - l)/2. Eq. ffTHI]) 
should be compared with Eq. fIST]) (with fi = (f) = 0), called now Eg, the energy calculated 
considering only a single state, \iPl) or \iPr), separately. For convenience we report here Eg, 
the energy of an unbalanced coherent-like state. 

Eg = ~NVl -z'^ + ^N{N - 1) (1 + z"^) . (151) 

Eq. (11511) has a single minimum in 2; = for U < Uc and two minima in z = Zo given 
in Eq. ( ]135p . exhibiting a sort of second order phase transition. The concavity at z = 0, 
(PEg/ dz'^\z=o = n{t + (A^ — l)f/)/2, in fact, is negative for U < Uc- The concavity of E at 
z = 0, instead, is always negative for any f/ < 0, (PE / dz'^\z=o = UN{N — l)/4. This means 
that there is always at least one finite minimum, i.e. Zo is finite for any U < 0. From Fig. [3] 
we observe a typical picture of a first order transition, with a jump in the global minimum 
at some U* < Uc, as shown in the left plot of Fig. HI For U > U* therefore, there is a 
minimum of the energy E given by 



while for t/ < U* the minimum is given by Eq. (I135p . In terms of Uc, Eq. (I130p . we have, 
therefore, 



2U 



o> u>u: 



z.^±{ V^^H^ ^ (153) 



However U* — )■ Uc for N ^ 00 (see Fig. |4]). More quantitatively, at leading order, we get 



At U = U* there is a small jump in the values of Zq, for larger than ten, which is given 

^■^ (1 ~ {V2N+i)^ y^^ ~ il^^^^y^'^ ^^'1 which goes to zero as {jfY^"^, in the large N limit. 

What we have learnt is that, only in the limit of — ?■ 00, which implies E Eg, the 
picture of spontaneous symmetry breaking is correct, while for finite N there is always a 
finite imbalance for any attractive interaction, but always with a null population imbalance. 

The number fluctuations on each site, ai^R = {n\ ^)(a — is given by 

N NiN-l)z'^ , , 

which goes from a = N/4 in the coherent-like state {z = 0) to a = in the NOON 



state {z = 1) |41|. 
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FIG. 3: (Color online) (First three plots) Energy E as a function of z, Eq. (|150p . for N = 300 and 
U = xUc, with X = 1.036,1.048,1.062 from left to right. There is a critical interaction U* < Uc 
at which the global minimum of E jumps from one value of z the another. The red dashed line 
is Eg as a function of z, Eq. (|15ip . for the same values of N and U. The energies from exact 
diagonalization, for those values of interaction, are, respectively, Eex = —227.97, —228.92, —230.05 
in units of t. (Last plot) Energy differences, between the global minimum of E and the exact 
ground state energy E^x, i-e. ^E^x = min(£') — E^x, (solid blue line) and between the global 
minimum of Es and E^x, i-e. ^E^x = mm(Es) — E^x, (red dashed line), for = 300. We observe 
that min(ii^) is much better than mm{Es) for U > U*, and deviates from exact ground level mainly 
for U ^ U*, although, the relative error at that point is still very small, less than 0.07%. 




-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 -1.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 -1.6 -1.4 -1.2 -1.0 -0.8 -0.6 



-U/U, -U/U, -U/Uc 

FIG. 4: (Color online) (Left plot) Imbalance Zo, calculated as the global minimum of E given by 
Eq. (I150|) . as a function of U in units of (—Uc), for N = 300 (solid line) and = 50 (dotted line). 
The red long-dashed line is the plot of Zo given by Eq. (jl35|) . i.e. the minimum of Es given by 
Eq. ()15ip . We can see that at U = U* < Uc a discontinuity of Zq occurs. (Central plot) Energy 
difference AEq = inm{Es) — mm(E), between the global minimum of Eg and the global minimum 
of E, as a function of U, for N = 300 (solid line) and = 50 (dotted line). AEo is always 
non-negative, large for > U > U* and almost zero for U <U*. (Right plot) Second derivative of 
the ground state energy E^x, obtained by exact diagonalization with respect to U, i.e. d'^E^x/ dU"^ , 
as a function of C/, for N = 300 (solid line) and = 50 (dotted line). The energies E(.x have been 
rescaled for a better comparison: E^x = -E'ea;/(294.55 1), for A^ = 300, and E^x = Eex/{5.68t), for 
N = 50. 



Let us now consider the visibility. Its analytical expression can be read out directly from 
Eq. (11441) . therefore 

, Vr3?+(v/T^)"-' ', 
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which is a = 1 in the coherent-hke state while a — )■ for z — )■ 1, namely going towards the 
NOON state. 

Finally we calculate the decay rate in the presence of a weak coupling to an external 
environment, as described before. From Eqs. fll44p - (11491) . we get, for F^, the following 
equation, written in terms of the previous quantities 



F^ = (y\cL - crP + F^JIchp + a^Re[(cL - cr)c^ 



(157) 



and, analogously, F5 = a\cL — cr]'^ + F^J|c/jp — aYRe[(cL — Cij)c^], where a is given by 
Eq. ffT55l) . a by Eq. ffT56l) and 



N N\ 2n 



(158) 



For z = {a = 1, a = N/A) we recover the coherent-like result, F^ = N/4, while for z = 1 
(a = 0, cr = N'^/4) we have F^J: = N/2. For large N, the asymptotic form of F!^ is 



Af-s>oo 4 



> ^ (1 + 



(159) 



while, at finite and large A^, F!: has a minimum at 2 ~ l.Q/yN, where T^/N ~ 0.11, see 



Fig. [5l However, the important point is that F^ does not scale faster then A^. In the presence 

0.5 [ 




FIG. 5: (Color online) The contribution F^^ to the decay rate, rescaled by A^, as a function of z, 
for N = 300 (solid blue line, with a narrow dip) and = 50 (solid red line, with a broad dip). 
The dashed line is the limit — )■ 00, Eq. (I159p . 

of a bath coupled differently with the densities located on the two sites (c^ 7^ cr), instead, 
the decay rate F^, for U < U*, scales as A^^ with the number of bosons, since the leading 
term is driven by the number fluctuations, F^ ^ (7\cl — cj?p, and a ~ A^^ for U < U*, as 
one can check from Eq. (I155p . evaluated at Zo, Eq. (ll53p . 

Before we conclude this section a brief explanation is in order. The discontinuity of Zo 
at U* cannot be directly measured since Zo is not an observable, being not equal to the 
population imbalance, which is always {rii) — (n^j) = 0. The small jump in Zo, however 
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could produce a small discontinuity in the visibility or in the variance of on-site number of 
bosons as functions of U. By exact diagonalization of H, instead, those quantities do not 
exhibit any discontinuities. However, the exact results are in a very good agreement with 
our analytical results away from U*, while, close to U*, those quantities, like the visibility, 
for instance, obtained by exact diagonalization, are smooth. Nevertheless, their derivatives 
at U* are large, although not infinite. The analytical expression for the energy is also in 
perfect agreement with the exact numerical result E^x, see the last plot in Fig. [3] and Fig. 
m It behaves as ~ — tN/2 and ~ UN'^/2 — tN for, respectively, very small and very 

large interactions. The point of crossover is identified by looking at the mimimum of the 
second derivative of E with respect to U (see the right plot of Fig. H]) which coincides with 
U*, the point where (min(ii^) — min{Es)) ~ (compare the right plot and the central plot of 
Fig. H]). We have checked this result also for different values of A^. Therefore, this value of 
interaction, in the presence of a finite number of bosons, can be promoted as the crossover 
point, which characterizes the loss of coherence and, as we will see in the next section, is also 
the point where the quantum entanglement reaches its maximum. Finally, for U < U*, as 
explained in the following paragraph, the cat state becomes extremely fragile under a small 
offset between the two on-site energies, i.e. /i 7^ 0. This fragility is also revealed, as we have 
seen before, by the decay rate in Eq. (11571) . when the two sites are coupled differently to a 
source of noise induced by an external bath, i.e. cl 7^ cr, since, in this case, the instability 
is driven by the large number fluctuations which goes like N"^ for U more negative than U*. 

A. Breakdown of the cat state 

The fragility of the Schrodinger cat state can be shown not only by analyzing the decay 
time = 1/r^, which goes like for U < U* and cl 7^ cr, but also by energetic 

arguments. As one can see from Fig. H] (central plot) the energy difference AEo between 
the minimum of E, calculated in the cat state, and the minimum of Eg calculated for a 
coherent-like state with finite imbalance z, is always positive and becomes almost zero for 
U < U*. Actually, the positions of the global minima of the two energies almost coincide, 
below that critical interaction, see Fig. [31 However, the energy obtained using the cat state 
is unaffected by a finite value of the chemical potential difference /i between the two sites, 
unlike the energy of a coherent state. The latter is modified by an additional term fiNz. 
Therefore, for large and U < U* an almost infinitesimal value of /i makes the cat state, 
superposition of two unbalanced coherent-hke states, energetically unfavourable with respect 
to a single coherent-like state. In this case, the number fiuctuations and the visibility become 
simply those in Eqs. fITB]) . (1771) . or, in terms of z, 



0" = — 




(160) 
(161) 



a = 
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Notice that a drops the A^^ dependence and becomes hnear in the number of bosons. Finally, 
let us consider the decay rate. If /i < the cat state collapses to and the decay rate 
r^, in terms of is given by 

= J {|cl - crWI - z^) + \ch\\l - zf + 2v/r^Re[(ci(l + z) - c^(l - ^))c*j} . 

(162) 

For one has to interchange cl with cr {cl ^ cr) and change the sign of z (2; — )■ —z). If 
/i > 0, the ground state is, instead, \iPr), therefore, the decay rates are the sames as before, 
providing that z ^ —z. 



X. ENTANGLEMENT FOR [/ < 



A. Entropy 

After tracing over one site, for instance, the right one, we get 
1 „ N\ 



P 



TrK(|^(a)(?/'@|) 



^LL ^ ^RR ^ ^LR _^ pRL^^ 



(163) 



where p^'^' = T^Tr/j (|'?/'s)(?/'^|), are diagonal matrices with finite elements given by 



.0 ^ LL _ RR 

Pi — Pi — Pn-e 



>\2e\t<\2{N- 



Pi"" = PN-i 



(164) 



(165) 



where + = 1. Eq. fll63p is therefore a diagonal matrix, = Su'Pi, with 

1 / „ „ /N- 



Pe 



2 l+(Re[e*e<]) 



N 



Pi + PN^i + 



Re 



N 



(166) 



As already said, at the ground state = 2n7c, therefore cos0 = cos N(j) = 1, and remember- 
ing that l^^'^P = (1 ± z)/2, we have simply 



2(1+ {VT^) 



N 



Pi + PN-e + r,N 



where p^, in terms of z, is given by = 2^ ^ ^ j (1 + zY{l — zY^ ^\ It is clear, also by 

looking at Fig. |6l that the profile of the reduced density matrix has a bimodal shape, peaked 
at the values N{1 + z)/2 and N{1 — z)/2, with an additional interference term. When z is 



(167) 



27 




FIG. 6: (Color online) Reduced density matrix for = 50 and z = 0.1 (yellow dotted line), z = 0.3 
(red dashed line), z = 0.8 (blue solid line). The lines are guides for the eye. 

small the two peaks merge together, while for z close to 1 the two peaks are far apart. 
For z = 0, in fact, we recover the coherent-like reduced density matrix 

P. = ^(^), (168) 
while for z = 1, namely in the NOON state, we have 

Pe = ^ {5i,o + kN) ■ (169) 

Now we can calculate the entanglement von Neumann entropy, whose expression is here 
reported, 

N 

^ = - J]p,log2P£. (170) 
For z = and for large N the reduced density matrix approaches a gaussian distribu- 



tion. Of ~ , exp 

' ^ y/nN/2 ^ 



(e-N/2f 



therefore, the asymptotic behavior for the entanglement 



N/2 

entropy in the coherent-like state is given by Eq. (11131) . namely 



So = l\og,[^N] . (171) 



For z = 1, in the NOON state, instead, since the reduce density matrix is given by Eq. (11691) . 
and then the entropy is simply 

^1 = 1. (172) 

For finite z, in the large N limit, we can approximate the reduced density matrix as p£ ~ 
|(P£ + PN-i)y sum of two gaussians, therefore, if z is large enough that the two 

gaussians are well separated, the entropy has the following asymptotic form 

S,c^^\og,{27ieNil-z^)) , (173) 



28 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



FIG. 7: (Color online) Entanglement entropy, S, in the cat state, as a function of z, for N = 300 
(upper solid line) and = 50 (lower solid line). The dashed lines are the asymptotic behaviors 
given by Eq. (jl73p . for N = 300 (upper dashed line) and = 50 (lower dashed line), which match 
the exact entropy for z > 2.6 /^/N. The straight dotted lines are the values of S for z = 0, given by 
Eq. p7ip . i.e. So, the entropy of the coherent state, for = 300 (upper dotted line) and = 50 
(lower dotted line). 

where, of course, z'^ < 1 — 2/TTeN in order not to go below the NOON entropy or even 
to diverge. Approaching 2; — )■ the two separated gaussians of the reduced density matrix, 
see Fig. [6l which gives the entropy in Eq. fll73p . start to see each other when the distance 
between the two peaks, which is Nz, becomes smaller than a certain number, n^j, of standard 
deviations, cr^ = a/ A^(1 — -2^)/4, of each gaussian. Namely, the entropy reaches its maximum 
and then starts to decrease approaching 2 = 0, when z fulfils 

Nzc^2agn„, (174) 

whose solution is given by 

\/ni + N 

From Fig. [TJ and further checks, in fact, we have found that the number of standard de- 
viations, consistent with such a description, is 71^^ ~ 2.6, therefore the maximum value for 
the entropy is reached at zs^^^ ~ 2.6/-\/iV, corresponding to having a distance between the 
two peaks of the reduced density matrix approximately equal to five sigmas. If z were given 
by Eq. fll35p . the entropy would reach its maximum at Us^^^ — Uc\J ^^-jp^ — Uc (^l + 
However, since Zo\u;-e < -^^Smax < Zo\u*+e, then the maximum entropy occurs at 

Us^.. = U: , (176) 

where [/* is given by Eq. f ll54p . For large A^, zs,^^^ goes to zero, and, as a result, the 
asymptotic value of the maximum entropy can be written simply substituting z with in 
Eq. (UTS]), getting 

Sma. ^ \ log2 (27reAr) = So + l, (177) 
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which does not depend on the interaction strength. In conclusion, the entropy grows up to 
1 above its value at f/ = 0, Eq. ( 11711) . within a short range of interaction close to the value 
f/*, Eq. fll76p . Strikingly, from Eq. fll73p . we have that the entropy exceeds its coherent 
value {S > So) for < z < \/3/2. In terms of the interaction, this means that 

S> So, for 2Uc<U <(}. (178) 



B. Fisher information 



From the definition of Fisher information, Eq. f ll24p . and after verifying that 

{nLnR)m- {nL)@{nR)@ = -a (179) 
where a is given by Eq. fll55p . we find the following remarkable exact result 

valid along all the attractive regime, from z = {U = 0, coherent-like state), where a = N/4 
and F = 1/N, to z = I (U -oo, NOON state), where a = N^/A and F = 1. Even if the 
system collapses into an unbalanced coherent-like state, Eq. f llSOp is still valid, providing 
that 0" is given by Eq. f ll60p . Surprisingly, Eq. ( I180p is the same relation as that obtained 
also for repulsive interaction, Eq. (I128p . 



XI. SUMMARY OF MAIN RESULTS WITH PLOTS 

In this section we show some plots of several quantities: energy (Fig. [S]), coherence 
visibility (Fig. |9]), Fisher information (Fig. [TOj) . and entanglement entropy (Fig. [TT]) . all as 
functions of the interaction strength, both in the attractive and repulsive regimes. We used 

= 50 in order to show features visible only for not too large number of particles. 



XII. CONCLUSIONS 



We have studied, in the repulsive regime, the interpolation from a spatial separated 
Fock state, ground state in the limit of strong repulsive interaction of the Bose-Hubbard 
Hamiltonian, to a delocalized coherent state, ground state for free bosons. We have used 
as many-body wavefunction a simplified permanent state which allows us to write all the 
physical quantities we are interested in in terms of polynomials of the single particle overlap 
parameter which can be determined variationally. Although not exact, our single particle 
approach provides a quite transparent and intuitive description of the crossover between the 
two exact limits. We have calculated the energy, the charge fluctuations, the decay rate, 
the coherence visibility, the entanglement entropy and the Fisher information with such a 
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FIG. 8: (Color online) Ground state energy Eq = min(£') as a function of U for N = 50, for 
attractive (left plot) and repulsive (right plot) interactions. Left plot (attractive interaction): the 
energy is described by Eq. (|150p at Zo given by Eq. (|153p . Right plot (repulsive interaction): 
the energy is described by Eq. (|103p . In both the plots, the dashed line is the coherent energy, 
Ec = —t/2 + UN{N — l)/4, showing that a finite imbalance lowers the ground state energy. The 
ground state energy obtained by exact diagonalization, E^x, coincides, on the scale of the plots, 
with the analytical result Eo (solid line), in both the regimes. 
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FIG. 9: (Color online) Visibility a as a function of U for N = 50, for attractive (left plot) and 
repulsive (right plot) interactions. Left plot (attractive interaction): the visibility as in Eq. (jl56p . 
calculated at Zo given by Eq. (jl53p . The dotted line is the limit — )• oo. Right plot (repulsive 
interaction): the visibility described by Eq. (jl04p . 



trial wavef unction. We have derived a non-perturbative, simple analytical expression for 
the ground state energy in the large N limit. We have shown that the charge fluctuation, 
null in the Fock limit while proportional to for zero interaction, in the presence of finite 
interaction scales as A^^^^ with the number of bosons. Moreover we have shown that the 
decay rate, in the presence of a weak coupling with an external environment, scales as N^^'^ 
in the intermediate case, i.e. finite U, and that it is proportional to the inverse of the 
Fisher information. In this regard, the overlap of single bosons can be considered as the 
key parameter which contains the quantum information. Moreover we have shown that, 
for U <^ N, the Fisher information F and the on-site number fluctuations expressed by 
the variance a are related by the equation F = 4a /N"^. In the attractive regime we have 
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FIG. 10: (Color online) Fisher information F as a function of U, for = 50, for attractive (left 
plot) and repulsive (right plot) interactions. The Fisher information is described by -F = Aa/N"^ 
where a is the variance of the number of particles. Left plot (attractive interaction): a is given 
by Eq. ()155p . calculated at Zq in Eq. (|153p . The dotted line is the limit N ^ oo. Right plot 
(repulsive interaction): F given by Eq. (jl23p . calculated from full correlators, Eqs. ([30]) . ([32|) . with 
parametrizations Eq. ([95|) . at w = oJq, Eq. ([99]) . The dashed line, which almost coincides with 
the solid one, is simply given by Eq. (|126p at uJq- 
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FIG. 11: (Color online) Entanglement entropy 5 as a function of U, for N = 50, for attractive 
(left plot) and repulsive (right plot) interactions. Left plot (attractive interaction): S is given 
by Eq. (1170p . with p£ obtained from Eq. ()167p . calculated at Zq in Eq. (1153p . The dashed line 
corresponds to So, Eq. (jl7ip . 5 exceeds So for 2Uc < C/ < 0. Right plot (repulsive interaction): 
entropy calculated through Eqs. ([W]), with n = k = N/2 and Eqs. ([Ml), ([95]), at ujo, 

Eq. (IMI). 

calculated again the energy, the variance of bosonic numbers a, the decay rate, the visibility, 
the Fisher information and the entanglement entropy by using a symmetric superposition of 
imbalanced coherent-like states, namely a Schrodinger cat state. We have shown that this 
state, in the absence of a local offset between the two on-site energies, has always a lower 
energy with respect to a single coherent-hke state which usually is supposed to describe well 
the low interaction regime. As a result, the symmetry breaking for finite number of bosons 
does not occur, and the corresponding critical interaction Uc, shifted toward more negative 
U, i.e. U* ^ Uc{l + -^), is actually a crossover point separating two different regimes, the 
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coherent regime and the incoherent one; in the latter regime the cat state becomes extremely 
fragile. In the presence of an external bath, if the environment is coupled differently to the 
two local densities the relaxation time is very short, inversely proportional to the number 
fluctuations which go like N"^ for U <U*. Remarkably, we have found that the equation for 
the Fisher information, derived in the repulsive regime, i.e. F = 4a /N"^, is exactly valid also 
for the whole range of negative interactions. Finally, we have calculated the entanglement 
entropy and found that, close to the crossover point, U*, the entropy reaches its maximum 
value, which, in the large N limit, is equal to the entropy at U = 0, increased by one (notice 
that 1 is the entropy in the NOON state). Another important and final result, valid for large 
N, is that the entanglement entropy exceeds its coherent value (at U = 0) for 2Uc < f/ < 0. 
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